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A straightforward analysis slightly extending the work of Kawakami 
and Nishizawa and of Gloge and Marcatili provides some insight about 
the extent to which the mode differential delay in a graded-index fiber 
can be minimized. We show that an ideal fiber (no mode mixing) with 
uniform mode excitation and loss and uniform material dispersion can 
theoretically have an rms pulse broadening due to mode differential 
delay as small as about 0.02LA 2 no/c. We suggest that further im- 
provement can result through recognition of differential mode loss and 
by accurate control of the (non-zero) rate of change of dispersion with 
fiber index. 



I. INTRODUCTION 

It is known through simple first-order analysis that differential delay 
between the propagational modes in multimode optical fibers can be 
greatly reduced by grading the optical index of the core so that the 
index 

n = no(l - AR 2 ), (1) 

where R is the fiber radius normalized to unity at the core-cladding 
boundary. It is also known that even less differential delay can be realized 
theoretically by slightly perturbing the gradient from this parabolic 
shape. 

We have taken some direct steps based on existing analyses to deter- 
mine how much further improvement might be realized if the optimum 
gradient could be realized in an "ideal" fiber, where geometry is invariant 
over its length (no mode mixing) and where material dispersion is in- 
variant with radius. The approach is very simple and will be so stated, 
but the algebra is tedious and what little has been included will be found 
in the appendix. 
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In their very nice analysis, Kawakami and Nishizawa 1 showed that 
an improvement in fiberguide impulse response could be obtained by 
perturbing the parabolic profile through the addition of a small 
fourth-order variation in the index gradient. They suggested that the 
minimum pulse width would be obtained when the fourth-order coeffi- 
cient, 5, lies between the values %, where all meridional modes are syn- 
chronous, and 1, where circular spiral modes are synchronous. Minimum 
total pulse width, t, in fact occurs when 5 = %, and minimum rms width, 
a, occurs when 8 - % Their expression for n can be written: 

n = n [l " 2A# 2 + 5(2 A) 2 /? 4 ] 1 / 2 . (2) 

It will be seen presently that when the index gradient is near optimum 
it takes some care to keep track of which propagating modes are the 
fastest and which are the slowest. The overall pulse width, t, is simply 
the difference in arrival time between the slowest and fastest modes. If 
we plot tvs 8 [found by solving (5) for all mode numbers at each 8 and 
taking the difference between the extremes], the curve is continuous but 
its derivative is not (Fig. 3 below). Personick 2 has pointed out that the 
rms pulse width, a, (the second moment of the received pulse) is more 
useful for fiberguide system analysis than is t. 



r2 = 



C t 2 p(t)dt r C tp(t)dt-i 
C p(t)dt t Cp(t)dt 



(3) 



where p(t) is the power arriving at a given point at time, t. Note that 
da/d8 is continuous, hence minimum a can be found by direct compu- 
tation. 

Gloge and Marcatili have shown in their analysis 3 that minimum total 
pulse width occurs when 

n = n [l - 2Aft 2 < 1 - A >] 1/2 . (4) 

For convenience, we introduce a multiplier, p, in the exponent of R in 
(4), namely, 

n = » [1 - 2Afl 2 < 1 -' A >] 1/2 . (40 

By definition, minimum total pulse width occurs for p = 1; minimum rms 
width, a, however, occurs when p = 1.2. 

These results are found by determining in each case the time of arrival 
of energy propagating in the p, v mode (radial, azimuthal mode number) 
with respect to the arrival time of the zero-order mode: 

t ( d ) = m^? [(1 - 35/2)(2p + v + l) 2 + fi/2(i» - l) 2 ] (5) 

2cM z 
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RELATIVE POWER, p(t') 




Fig. 1 — Normalized time of arrival vs mode number (right) and mode count vs time of 
arrival (left) for n = n [l - 2AR 2 + (%)(2A) 2 fl 4 ] 1/2 ; i.e., 5 = %. 



t(p) = 



LA 2 n u 
2cM 2 



[(2/i + v + l) 2 - P M(2n + v + 1)], 



(6) 



where L is the length of the fiber, c is free-space light velocity, and M 
is the largest guided-mode number. 



In both cases, 

M = (2n + i> + l) max 



irno(a/\ )V2A , 



(7) 



where a is the core radius and Ao is the free-space wavelength of light. 
Also, since we must include both E and H waves and two polarizations, 
a fiber with near-parabolic core gradient carries a total of about M 2 
modes. 

Equation (3) can be solved easily if we assume all modes are equally 
excited [p(n, v) = 1] and integrate over all modes. So 

J-»M r {M-v)l2 r M nKM-vSH , 

t 2 {n,v)dpdv r 1 I tip, v)dndw{ 

o Jo Jo Jo I 

(J<- = 



r 2 = 



M 2 /4 



M 2 /4 



(8) 
Substitution of (5) and (6) into (8) (assuming M » 1 to simplify the 

algebra) and minimization with respect to 8 and p, respectively, produce 

the results already stated. 

We can substitute these minimum values back into (5) and (6) and 

plot arrival time as a function of n and i>, as shown in Figs. 1 and 2. Also 
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shown in these figures (on the left) is a plot of relative power vs time of 
arrival, again assuming uniform excitation of the modes. 

Figures 3 and 4 show t' and a' vs h and p, respectively. The prime de- 
notes normalization with respect to L A 2 mq/c. 



RELATIVE POWER. p(f 




Fie. 2 — Normalized time of arrival vs mode number (right) and mode count vs time of 
rival (left) for n = n [l - 2A/?2< 1 - 6 ^ 5 >] 1 / 2 ; i.e., p = %. 



arrival 




Fig. 3 — Normalized total pulse width, /, and rms pulse width & vs h for n = no[l - 2AJ? 2 
+ 6(2A) 2 fl 4 ] 1 / 2 . 

722 THE BELL SYSTEM TECHNICAL JOURNAL, MAY-JUNE 1977 




Fig. 4 — Normalized total pulse width, r', and rms pulse width & vs p for n = nnll — 
2A(R)« l -" A >] 1/2 . 

It is natural to question whether we can improve on these independent 
functional perturbations from the parabola by an appropriate combi- 
nation thereof. Combination of (5) and (6) yields: 



«<>-i£?[('-!) «>*-*»' 



] 



+ -(*>- l) 2 - pM(2 M + * + 1) |. (9) 

Substituting (9) into (8), as before, and minimizing with respect to both 
8 and p produces: 8 = %, p = %. Substituting these into (9) and neglecting 
the Is (large mode count) yields: 



t 



LA 2 M 
cM 2 



[m 2 + W + Q v 2 - Q M(2 M + r)l. (10) 



This is plotted in Fig. 5. 

Now combine (2) and (4') and let 8 = % and p = % to find 

n = n [l - 2Afl 2 < 1 - 2A / 3 > + (%)(2A) 2 £ 4 ] 1/2 



(11) 



as the near-optimum index gradient. Equation (11) can also be writ- 
ten 

n « n [l - Aft 2 - A 9 c(A)], (12) 

which may be more convenient if we are seeking the excursion of the 
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t" = -0.08 CONTOUR 



RELATIVE POWER. p(f) 



- 




o'- 0.0216 



Fig. 5— Normalized time of arrival vs mode number (right), and mode count vs time 
of arrival (left), for n = n [l - 2Afl 2 < 1 " 2 */ 3 > + (l/3)(2A) 2 fl 4 ] 1/2 ; i.e., p - %,« - %. 



improved gradient from the inverse parabola of (1): 

*(/?)« (^)fl 2 In (l/fl 2 )-y. 
This is plotted in Fig. 6. 



(13) 



II CONCLUSION 

The rms pulse width in an ideal fiber with parabolic gradient is 

LA 2 rto 

o-« 0.145 -. 

c 

For the profile determined in this paper, 

L A 2 /io 
c 



a = 0.0216 



(14) 



(15) 



Hence, nearly an order-of-magnitude improvement in fiberguide impulse 
response over that produced by an inverse parabolic gradient could be 
realizable in an ideal multimode fiber. 

Uniform excitation of the modes was assumed for this analysis. If the 
low-order modes tend to carry more power, 2 even smaller a could be 
realized and Fig. 1 would suggest that p -+ and 8 -*■ % might be more 
nearly optimum. 

Olshansky and Keck 4 have shown that when material dispersion varies 
with radius, the optimum gradient is significantly different from para- 
bolic and from those discussed here. On the other hand, it can be shown 
that procedures similar to those indicated here can be carried out; and 
if the dependence of the dispersion on optical index could be controlled 
to a predetermined practical value, even less differential delay than that 
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Fig. 6 — Index perturbation from inverse parabola for minimum rms pulse width [see 
eqs. (12), (13)]. 



indicated here could be realized. This does not necessarily contradict 
the results of the calculation by Arnaud and Fleming 5 based on an 
analysis 6 which includes tunneling modes at equal weight with the more 
clearly guided modes. They showed that if we assume a very specific 
non-zero dependence of dispersion on index, namely, that resulting from 
the inclusion of germania in silica, the differential delay of the optimum 
graded-index fiber is considerably degraded. We suggest, however, that 
if we carefully choose materials to enhance the index (a propitious mix 
of germania and phosphorous oxide, for example), we might provide a 
dispersion vs index dependence that would improve rather than degrade 
the mode differential delay, at least at a particular light wavelength. 

This goes well beyond our present ability to make measurements and 
control materials, however, and only leads us to conclude that im- 
provement in technology can potentially bring significant improvements 
in the information-carrying capacity of low-loss graded-index fibers. 



GRADED-INDEX FIBERS 725 



APPENDIX 

It is our purpose first to justify (5) and (6). Equation (5) comes from 
Ref. 1, eq. (29), which translates directly to 

c 2 c a 

9A 

- 5[6m 2 + M* +D + ("+ D(i» + 2)] — , (16) 

a 1 

where a is the core radius. 
This can be solved for /? to find 

c a 

cA 



a>n a z L\ 2/ 2 J 

The largest propagating mode numbers can be found by recognizing 
that the phase constant, 0, can be no lower than the phase constant, k a , 
of an unbound wave in the cladding, where 

n a ~no[l-2A]i/ 2 ~n (l-A). (18) 

So 

/3 min «^-(2M + ,+ l)^^ = ^«^(l-A). (19) 

c a c 

Eq. (19) yields 

(2, + , + l) max = M = ^^ (20) 

M2 = ^^. (21) 

Substitution of (20) into (17) shows that each term on the right-hand 
side thereof is of order A smaller than the previous. 

The time of arrival of energy carried in any mode through length, L, 
of fiber relative to the arrival of energy in the zero-order mode is 

tmL \m_m ( 22) 

L da) dco J 

Differentiation of (17) and substitution of that [utilizing (20)] into (22) 
produces eq. (5). 

Equation (6) is found from substitution of GM(17) [i.e., eq. (17) of 
Gloge and Marcatilli in Ref. 2] into GM(19), recognizing that the a of 
Gloge and Marcatili is our 2(1 - pb). First, we must identify m/M of 
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Fig. 7 — All modes that lie under the line intercepting the n axis at M/2, the v axis at M, 
are propagating modes. 

GM(17). The symbols are unfortunate, but we find this equates to our 
(2/x + v + 1) 2 /M 2 . 

This is easier to justify than derive. All modes are associated with 
positive mode numbers that lie in the n, v plane, as shown in Fig. 7. We 
can plot the identity of (20) on the plane to bound the crosshatched area 
in Fig. 7 within which lie all propagating mode numbers. The total 
number of propagating modes is four times the total mode-number 
combinations in that area, which can be seen by inspection to be M 2 . 
Gloge and Marcatili 2 showed that for their assumed functional index 
variation, the phase constant depends on the total mode number inde- 
pendent of the ratio of n to v. If, then, we draw a dotted line parallel to 
the limit line, as in Fig. 7, the sum of all modes with phase constant 
greater than that represented by those on the line are identified by 
numbers that lie within the double crosshatched area. This Gloge and 
Marcatili call m, and we call (2/j + v + l) 2 . Their M is the total number 
of propagating modes, which is our M 2 . [We have neglected the extra (1) 
since we have assumed M » 1. We would have to be more precise if only 
a small number of modes were involved.] 

We can now translate GM(17) to be 

(2m + * + 1) 2 ,» /M2 
or 

«o M - A(8 »t; +1> - < 23 > 

M 
This in GM(19) results in (6). 
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Performing the integrations for (8) is simple but tedious. A helpful 
formula, however, is 

J*(M-,)/2 jtfn+l-^n+l 

o 2(n + 1) 



The evaluations of (8) for the two cases of (5) and (6) are 

LA 2 n 1 fl 25 , 75 2 ~|i/2 



(25) 



For the combined case, 

Minimizing these is straightforward. 
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